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ABSTRACT 


A method is derived for finding the linear response and 
loading transfer functions for the lateral aerodynamic case 
of airship flight through atmospheric turbulence. The func- 
tions obtained are in a form that can be applied to the 
various spectral analysis methods used to predict surviva- 
bility currently employed by designers. A numerical example 
using the USS AKRON (ZR-4) is presented. The results show 
that peak motion response and loading occur when the en- 
countered spectral component has a wavelength equal to the 
airship length, and that simple feedback of heading angle 


does not significantly decrease this peak. 
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i ene ODUCT ION 


When a close look is taken of the history of airship 
flight, it becomes evident that turbulence is a prime cause 
for concern in the design of lighter-than-air (LTA) craft. 
The spectacular crashes of the airships Shennandoah, Akron, 
and Macon are perhaps the most obvious reminders of the power- 
ful effect the wind can have on these fragile vehicles. [It 
is imperative that proper consideration be given to the gust 
response of an airship, both dynamically and structurally. 

At the time the great rigids were built, designers had 
only a cursory knowledge of turbulence and how to deal with 
it. Burgess [Ref. 1], for example, dedicates only one page 
to the subject, and that is mainly a warning that gusts en- 
countered in flight can place larger loads on an airship 
than any maneuver of which it is capable. He suggests using 
the standard gust analysis technique of the day--a "fixed- 
in-space" flight through a ramp shaped gust. This was the 
method used to evaluate all aircraft at the time, and pro- 
vided some measure of confidence that the structure would 
withstand the stresses of flight. Its principal weakness 
was that it took into account only the direct gust loads. 

Since that time, the methods used to analyze response to 
turbulence have been greatly improved. Over time, the ramp 


shaped gust gave way to other gust shapes, finally settling 


bs 





on the (l-cosine) gust as standard [Ref. 2]. By the mid- 
fifties enough data had been gathered on the nature of the 
turbulent wind to derive some basic statistical relations. 
This, in turn, allowed the development of the power spectral 
method for aircraft load analysis (see, for example, Press 
and Meadows [Ref. 3]). The spectral method has been applied 
to many types of aircraft over the last twenty years. At 
present, both methods are used in analysis of flight struc- 
tures, the one giving the more conservative structure de- 
termining the final design. 

While the development of turbulence modeling techniques 
progressed, applications to airship technology were slower 
in coming. Because the success of heavier-than-air (HTA) 
craft made the slower airship economically less attractive, 
LTA research lagged. In the period between World War II and 
the Arab Oil Embargo of 1973, the only significant contribu- 
tion to the study of airships in turbulence was by Calligeros 
and McDavitt [Ref. 4]. This paper presented a method of 
analysis that allowed a stable airship to respond dynamically 
to both sinusoidal and (l-cosine) gusts. Thus, the inertial 
and aerodynamic reaction forces of gust encounter could be 
included in the model. Further, by using the sinusoidal 
representation for gust shape, it was possible to avoply spec- 
tral methods to the analysis. 

The recent interest in LTA brought about by rising fuel 


costs has increased the research in all aspects of airship 


16 





Flight. The spectre of the great rigids disintegrating in 
turbulence makes it imperative that the designer have ade- 
quate means to predict an airship's response to gust pene- 
tration. Current research is aimed at supplying that means. 
DeLaurier and Hui [Ref. 5] refined the technique of Calli- 
geros and McDavitt [Ref. 4], by introducing refinements to 
the aerodynamic cross-flow model [Ref. 6] and allowing for 
emeerlity augmentation through pilot control input. The 
model allows statistical prediction of an airship's dynamic 
response and operational lifetime for various combinations 
of speed, altitude and control gain. 

DeLaurier and Hui's work (as well as most others dealing 
with this subject) concentrated only on the longitudinal 
aerodynamic case. This thesis proposes to apply their model 
to the lateral case, enabling the response to side-force to 
be calculated. Bending and shear in the horizontal plane, 
as well as twisting moment, can then be taken into account 


when predicting airship life expectancies. 
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i. Lait TURBULENT WIND 


The motion of the atmosphere is very complex. Shearing 
stress between layers of different speeds and at the ground, 
thermals caused by solar heating, weather fronts, vortex 
shedding behind obstructions and aircraft, plus many other 
phenomena, all contribute to a velocity field that is most 
difficult to describe. Of the methods available, that chosen 
will depend on the purpose for which it is used. A goal in 
the design of any flight vehicle is safety, with performance 
adequate for the mission. This dictates using models giving 
reasonable estimates for the design parameters, and it may 
Or may not be necessary to closely match the physical reality 
to do this. In any case, the designer must be familiar with 
the advantages and limitations of methods available in order 
to choose wisely. What follows is a brief review of some of 
the turbulence models currently in use and how they apply to 


airship analysis. 


pee, LHe DISCRETE GUST 

As mentioned in the introduction, the discrete gust model 
has been used to analyze aircraft for many years. It is 
especially good when response to the passage through a steady 
velocity gradient, such as a thermal, mountain updraft, or 


jet stream, is desired. The method has been improved 
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steadily until, as Etkin [Ref. 7] points out, it has attained 
a high degree of sophistication. 

Figure 1 shows the shape for a (l-cosine) gust, where ae 
and dq. are maximum gust velocity and distance along the 
flight path of this maximum. By varying these parameters, 
the gust severity can be controlled. The value for do is 


prescribed by the Federal Aviation Administration (FAA) as 
2d_ = 25¢ 
m 


The size of a 1s dependent on airspeed and altitude, and is 
shown in figure 2 for three values of equivalent airspeed. 
The 25c wavelength was chosen because it historically couples 
with the short period pitch mode of a rigid aircraft to pro- 
duce the largest load factors. Calligeros and McDavitt [Ref. 
4] showed that, for airships, the maximum loads occur when 
the wavelength is equal to the airship length. 

tie Bertvsh dictate (ARB CAR CH D3-3) that the gust para- 
meters be chosen to produce the peak response with aircraft 
flexibility taken into account. In this way, the model is 


"tuned" to the aircraft, thus assuring a conservative design. 


B. RANDOM TURBULENCE 

Extensive measurements of the atmospheric velocity field 
have been made, and the techniques involved are well estab- 
lished and reliable [Ref. 2]. They show that the velocity 


vector is best characterized by a random function of space 


I, 





and time that is, in general, non-homogeneous, non- 
stationary, and anisotropic. The exact function has not 
been developed due to its tremendous complexity. Until 
enough data is collected (1f ever) to allow the precise 
formulation, certain simplifying assumptions must be made 
to enable flight vehicle analysis. 

One assumption that applies everywhere except in the 
planetary boundary layer (below about 1000 ft), is that the 
turbulence is 'homogeneous', that is, the statistics of the 
field do not vary through space. In the boundary layer, 
scale length and intensity are homogeneous in the horizontal 
plane, but not vertically. Another assumption made is that 
the turbulence is 'stationary', or statistically time con- 
stant. Over the time periods of interest to flight this 
approximation is quite adequate. Also, the turbulence is 
assumed to be isotropic (again, except in the boundary 
layer), making the statistics invariant with orientation. 

One last simplification used to model atmospheric tur- 
bulence is the ‘frozen field' or Taylor's hypothesis [Ref. 
8]. The change in the velocity field perceived by an air- 
craft as it passes with speed Us through the air 1s given by 


the. substantial derivative 











Dee): .. 2) o( ) 
I se eo ime 
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Taylor's hypothesis states that for all but the smallest 
values of U5: the second term dominates and the first may 
be ignored. The result is that the correlations and spectra 
reduce to three-dimensional functions of space only. Phy- 
Sically, this means the velocity field is 'frozen' in time, 
and the changes are due only to displacement. 

In an effort to specify an acceptable lower limit on U_, 
for which the frozen field applies, Dobrolenskiy [Ref. 9] 
cites studies comparing records of turbulence spectra gath- 
ered by captive balloon and an aircraft flying nearby at the 
same time. Within the margin of error, the two are statis- 
tically quite similar, and he concludes the lower limit on 
Us 1s comparable to the convection velocity (for all prac- 
tical purposes, the mean wind speed). Etkin [Ref. 7] points 
out that the vehicle speed can be as low as one-third the 
wind speed for good results. Note that the only vehicles 
capable of less than this velocity are LTA and VTOL craft, 
and then only when they are convected downwind with the air- 
mass. In hover, or upwind flight, the hypothesis holds. 
For that small portion of the flight envelope in which it 
does not, the forces generated on an airship's structure are 
small and, therefore, do not present a problem. 

The techniques for dealing with isotropic, frozen turbu- 
lence are well known [Ref. 8]. For those not familiar with 


the mathematical background necessary to deal with the 


a 





Peesece in Geptn, Chapters 2, 3, and 13 of Etkin [Ref. 10] 
will provide a good primer. 
ime Turbulence Organization 

Anyone who has seen leaves swirl on an autumn day, 
or watched someone blow smoke rings, has an intuitive under- 
standing that the motion of the atmosphere is not completely 
arbitrary. What happens at one location effects conditions 
at another. Fuidynamicists characterize this interdependence 
by using expressions relating the stress and strain in the 
fluid. The technique has been applied to turbulence [Ref. 
11] with some success in predicting the actual velocity field 
Of a boundary layer type flow. 

For flight vehicle analysis, a more convenient method 
of specifying the velocity field is the spectral decomposi- 
tion of the three-dimensional homogeneous vector field [Ref. 
12}. Figure 3 shows an aircraft flying through a (two- 
dimensional) sinusoidal wave of shearing motion. The velo- 
city change from the mean is given (for the lateral gust 
component) by 


L(Q25x' + Qoy") 
ay (xi, y') =e dc. (aL) 


where 24 and 2 


y' directions respectively, and C5 1s the complex amplitude 


are the wave number components in the x' and 


of the lateral component. If the vehicle penetrates the 
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erela with velocity U5: the coordinates become x' = x + Ut. 
and y' = y in the body fixed system. Equation (1) then 


becomes 


dv (X,Y) = 2 le es ——— dc, (2) 
The air velocity over the vehicle is then periodic with wave- 
lengths (27/2) 5) and frequencies (27 9U,/2). The total 
field is made up of the superposition of these spectral 
components, much as a Fourier series represents a random 
Sealar. 
Pea woon Dey Pistribution and Spectra 

With the expression for a single spectral component 
available, the next difficulty is determining the probability 
distribution of the individual frequencies. The power spec- 
tral density of a time varying function, X(t), is defined 


(in terms of wave number) as 


AL 


6(2) = lim ao | Bite, oe) neice (3) 
AQ+O 0 
Too 


where $(2) is expressed in (ft/sec) */(radians/ft) , and T is 
the duration over which X(t) 1s measured. The value is usu- 
ally computed by taking the autocorrelation function R(T) 
(Ref. 10: Chapter 2], and performing a Fourier transforma- 


tion, thus 


re) 





ote = s | R(t)er Taz (4) 


Figure 4 [Ref. 2] shows the spectra of three different 
meteorological conditions. While varying in detail, they 
all exhibit the same decreasing trend at higher frequency. 
The vertical dimension is a measure of the intensity of 

the turbulence at that particular frequency, and the square 
root of the area under the curve a measure of the overall 
rms gust velocity [{Ref. 13]. As a point of practicality, 
only the area under the actual measured curve is included, 
because values of the spectrum at higher frequencies con- 
tribute little to the response of an aircraft. 

From analyzing large numbers of samples, it iS ap- 
parent the probability distribution is non-Gaussian [Ref. 
14], with very high and very low values of intensity occur- 
ring more frequently than predicted by a normal distribution. 
However, the vast majority of values do fall on a Gaussian 
curve, so it is reasonable to use this assumption for most 
applications. This is very beneficial because, whereas 
Gaussian input to a linear system results in Gaussian out- 
put, response to non-Gaussian input is, in general, unknown. 
In the analysis of flight vehicles, linear system models are 
used extensively, making the assumption of normal distribu- 
tion most desirable. [In order to account for the large 


gusts omitted by this type model, the (l-cosine) method can 
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be employed. This is the practice recommended by many cer- 
tifying agencies. 
Two Gaussian models in current use are the Dryden 


spectrum 





‘ ieee, + 30-0° 
ee te nn 22,2, 2 (5) 
ees 9) 
and the von Karman spectrum 
Pee 1 339 Ra)i- 
cama): ae : (6) 
33 7 aw AN 


fl + (1.3390) ~) 


The first was developed to define turbulence spectra in wind 
tunnels. It 1s the simpler of the two, but not as accurate 
as the second. For that reason it is not used as much today 
as in the past, but is mentioned here due to its historical 
Significance and the large number of references to it in the 
literature. 

Today, the von Karman spectrum is used almost uni- 
versally. In equation (6), oc is the rms turbulence intensity, 
and L is the "scale length"--a measure of the average eddy 
Size encountered. Testing has shown that the model is a 
reasonable fit to all levels of Eanes | Figure 5 isa 
plot of one set of experimental data along with the predicted 
turbulence spectra for severe storm conditions using various 


values for L. As can be seen, the agreement is quite good. 
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The values for o and L are variable and. appear to be 
functions of altitude. In addition, two standard categories 
of intensity are defined--"storm" and "non-storm". Table I 
(Ref. 15] lists values for non-storm (b,) and storm (b,) 
intensities, as well as scale length, currently used by NASA 
for horizontal atmospheric flight. In this table, the values 
Py and P> are the probabilities of encountering non-storm 
and storm turbulence, respectively, at the altitude speci- 
fied. Note that the values of : Given £Or Llight below 1000 
ft are representative values that are probably low, giving 
conservative (high) numbers of load exceedances per unit 
length of flight. These values are usable for structural 
analysis, but inappropriate for control system studies of 
Flight simulation where the vertical inhomogeneity must be 
taken into account. 

Equation (6) is the expression for the transverse 
spectrum needed to analyze the longitudinal response to 
vertical gusts. When dealing with lateral aerodynamics, such 
as the analysis done in Chapter III, it 1s necessary to use 


the longitudinal spectrum 





3°L 2 
© (Q) = ~ (7) 
tl 1 fl + (1.3392Q)2}>2/° 


where o and iE are the same as in the transverse case. The 
term “longitudinal" refers to the variation in gust velocity 


parallel to the direction of the mean wind, whereas 


26 





"transverse" refers to the perpendicular direction (both 


vertically and horizontally), thus 


ee 33 


It should be obvious from this discussion that flight direc- 


tion relative to the mean wind becomes important when dealing 


with lateral aerodynamics. 


| 





Pelee OF ANALYSIS 


ee aes MODEL 

The aerodynamic model used in this analysis is the same 
as used by DeLaurier and Hui [Ref. 5], except that it is ap- 
plied for the lateral case. The assumptions are as follows: 


1) the vehicle is perfectly rigid, flying at a reference 
velocity US through a constant density o 


ani) the motions are described by the lateral case only 
iii) the control provided by rudder deflection is linearly 


proportional to the yaw angle (y), that is, 
ve = KAY 


we 
This last assumption is the case for a helmsman using greater 
control the farther off heading he is perturbed, an assump- 
tion in keeping with operational practice. 

The turbulence model used is the von Karman spectra 
described in Chapter II. It 1S assumed that the turbulence 
1s composed of horizontal gusts only--either ue Or Ve de- 
pending on airship heading relative to the mean wind direc- 
EON. a 1s shown in the analysis.) 

i )eOue@en fxom Turbulence Components 

From Jones and DeLaurier [Ref. 6], the normal force 
on a hull segment of length d& is (see figure 6) 


dA 


eel: 2 8 
Ea ZOU, (Kees aene 2.0: ) cos| 5 ae 


5 5 cae: (Cae). g@in@ sin|9; 2rdé] 


(8) 





where: K is the hull potential cross-flow factor [(Ref. 6] 
@ is the angle between the hull centerline and Oe 
ris the radius of the hull segment. 


Differentiating with respect to §& to obtain a perturbation 


equation gives 


9 


Z 2 , Rees (9 A 
aF, = ZeU, | $ sin (26) sin|3 ~ 2 K cos(28) cos || ae 





| 
(Cq), (cosésin| 6| + sinde@s'| 6 | ad 2r dé | 


By limiting the analysis to the lateral aerodynamic case 
only, (9) becomes (8), the sideslip angle, and Fy becomes 
Yue the hull sideforce. Further, if (Bo). the undisturbed 
value of sideslip, is assumed to be zero--the usual case-- 


the above expression becomes 


R t 
Real: 2., dA O = 
(azn = ZOU, de 2008 (23,)cos =| lee, d G 
_ 2a Ge 
— pu, K ae aé ae 


Finally, for smail values of er 


“g 
dg == 
O 
and 
2 adie. oc 
(ay oh = 0U, K ag & Ty (9) 


oo 





The stabilizer forces are given by 


2 
U (v_) 
_ Oo S 
ee ee ene atk.) sh (10) 


where: H(k.) is the generalized Sears function as given 
by Filotas [Ref. 16] 


k= , the “reduced frequency" of the fin 





a) 1s the gust velocity at the fin mid-chord 


The propellers used to drive the airship vroduce a 
Side force when acted by the turbulence [Ref. 17]. Each 
thruster contribution adds to the total force and moment 
produced, and can be described, for the jth thruster-rotor 


Sembination by the following: 


2 
(Nee)... 
7 oO oe 
Wo)p, = 27 Sry (Cy) a, (11) 
O 
hg) p, = Mg) os Pom = nq) (12) 


Bem@etions (12) and (13) assume that the rotors are arranged 
symmetrically about the x-z plane, so that moments due to 


rotor offset in the y-direction cancel out. 
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2. Aerodynamic Forces and Moments Due to Airship Motion 


eae any- vic) 
en | US Kae u. 





One) : 
~ onl, se + Ugtk, | ae (14) 
where Koy Ky are the horizontal and longitudinal apparent 


mass coefficients respectively and 


Vie) = = r{ (la, = 


) 


Jt 


For the fins we have 





: 2 , wie 8° cy cv. 7 
(Y).=-p-o S. C ee | + [c | (HES) 
s — s| Ve ene | | Yaieecu YB). 20.7 | 
where v. =v- r (lan ~ ee 
wees YH a - 1.) 
(Ls = ie (i (16) 
ac — 
=e Z —- f ar 
Bees 29% 0 3s C51 Cnr}. 2U ie 


The superscript ac indicates the quantity in parentheses is 
taken about the fin aerodynamic center. 


Thruster forces and moments are given by: 


; ue (Vp) 5 
wt a ‘“yg) U (18) 
O 
alg = Yang (Bog ~ bes) (19) 


Ju 





Nw. = (Las Com ~ bay) (20) 


where 


(vp) =v- r(1 on as =a a DU a = one) 


3. Inertial Reaction of Airship EO se rodynamic POuces 


and Moments 

The "forces" covered here are those arising as reac- 
tions to airship motion. They are, in general, the negative 
of the forces and moments that give rise to the indicated 
airship translational and angular velocities and accelera- 
tions so as to produce a state of dynamic equilibrium 
(F - ma = 0). 


For the hull: 


(cree = -y' (6) (d_), (21) 
where y! =v + “<n - €) + Ur 

(aL Jy = =e os + ee (22) 

(dN Jy = -diI_4t + div5P (23) 
ror the empennage: 

(ts = “Mile “Ss (24) 

(Nis = -(I0) ot u (Io) gP ) 

(Ls = Ae hae , (ee eo) 


Se 





where dil! ele and di. are the moments and product of 
inertia respectively of the differential element under con- 
sideration, including all structure, air and gas contained 
in the airship. 


4, Bouyancy and Control Terms 


Referring to figure 7, the force due to bouyancy is 


given by: 
(Yo) 1 = ~(gdm ~- pgAdé) sind cosa. 
6 = roll angle 
2 = steady state angle of attack 


Differentiating to obtain a perturbation equation gives 


(dy, = (o9gAd& - gdm) cosé. cosa. do 
and letting %% = 0 results in 


(dye = [ogAdé - gdm] cosa, do (27) 


Control force is assumed to come from rudder deflec- 


tion, and acts through the fin aerodynamic center. 


Seve 7 Ko 
2 
O 


oa apm cane (28) 


5. Shear Force, Bending and Twisting Moment 


The hull's shear force at station (1) is obtained by 


summing the sideforce values from the nose, up to (1): 
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1 
Sel) | (€¥) + 
0 J 


iro 


. 
“%q) 9, + WY) pe (29) 


L 2 


and (a) is the number of rotors forward of station (1). 


Likewise, the bending moment at (1), measured along 


the centerline, is 


ile 1 a 
BM(1) = | eee, + | (dN )n + 2, Geta) (gh pt wg] 
0 0 J= : 
(30) 
Finally, the twisting moment at station (1) is 
1 1 a 
TM(1) = -| fen) (a2) | (dim), + | (Alug) hn 
0 0 0 
a 
- hm.) | (Y + (Y 
oe t,) (i gt, + wd (31) 


The term (dog) h is the torque contribution due to the center 
of gravity being offset from the central axis (see figure 7). 


teeas calculated from 


(dL dh = Dom dm 4 cosa, 


— 


B. FLIGHT DYNAMICS 


1. Bymamic Stability 


The equations for Lateral Dynamic Stability are taken 


from DeLaurier et. al. [Ref. 18] and given below: 


ACy sero + ae - (B - mg) cosa. = 2u(D8 + 2X) (32) 
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ACnsero + Song ~ *,B cosa,? = I,,Dr - I, Dp (ers ) 
Siero TEL. Ft 2B cosa,d = I,,Dp - I, Dr (34) 
Dd =H +f tana. (25) 
Dy =f seca, (36) 
Tear = Cy 28 - Cy3D8 + Cyt + Cy ,Dr + yoo + CyaUP (37) 
ACnaero = Cng? + CngD8 + Cn, % + CnzDF + CnpB + Cn5DB (38) 
Ps. = C128 + C13D8 + C1,F + Ci ,Dr + C1 yP + Ci5PP (39) 
Gy. & Kav (40) 
[l. - l__] lee alae) 
dCn, = = ie - = Sa ko (41) 
e e 
ats I coo ay gee 5) 
eee ee Ac. = ae kk (42) 
lis Sate Yo ee, Cc 
Cc e 


In this analysis 8 = a. = v7, thus equations (32) to 
(39) become . 
ACyserg + SCy, - (B - tig) cosa, = 2u(D + 2) (43) 
ieee SC, — 2,8 coso.$ = 1 (pt = 1 bp (44) 
SClaero + AC1e - 2,8 cosa, = I,,DB - 1D? (45) 
Dd=p +f tana, (46) 
Di: a F seca, (47) 
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AG 


eo = Cy av = Cy3Dv > Cy? - Cy,DF 


ee, = Cn 2¥ ct Cy DV + Cy f + Cn dt 


i _ = Sia’ = Cigbv 7 Cit + Cy, .br 


These equations are linear, and along 
through (42) can be written in matrix 


—~ 
V 


F¥> 


[D] 


) 
~K> 
i 


o> 


=> 


| 
| 
7 
where: v is translation in the y-dir 
Eis the rotation rate about 
me isethne rotation rate about 
¢ is the roll angle 
y is the yaw angle 


The matrices [C] and [BD] are given by 


[(C] = | Cy, = (Cy, -2u) Cy, tig 
Che Cn, Cny = 
Cle Cir 7 Cly ~ 2) 

0 tana, 1 

i 0 Seca , 0 
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+ Cy? + Cy5DP (48) 
+ Cepoe Co <Dp 49 

NpP npvP (49) 
+ ip? + C15PP (50) 


with equations (40) 


form as: 


ection 
the z-axis 


the x-axis 


-B) cosa, =e / 
. ae] ae) 
B cosa vs Tom’, | 
O c= Cr} 
“A +h | 
B cosa cm | 
O — k 
c cS 
0 0 
0 0 





[D] = ya! “Cys See 0 0 ; 
nis (Io “n2) = (I tCny) 0 0 

i 

“C15 lee (Ty f19) 0 a 

0 0 0 1 0 oH 

| 

ie 0 0 0 0 1 | 


To solve this system of equations, assume a solution of the 


form 
eS agl S* So =.Be oe — Re tot 
(22) 
aero are ) we rot 
which, when substituted into equation (51) becomes 
V 
R 
[Cc - ioD] | B | = 0 (53) 
| 9 
| 
io 


where ¢ is the non-dimensional stability root. This is an 
eigenvalue problem. A computer solution is performed to find 
the eigenvalues (stability roots) and eigenvectors (model 
vectors) for the control gains considered. For the subse- 
quent load-response analysis, dynamically stable cases must 


be chosen. 


3) 





2 eluepwLence Forcing Functions 


As mentioned above, the model chosen for the turbu- 


iemee is sinusoidal with Gaussian statistics, in particular, 





Ven) cosa_)} 
+ = [ exp|iwt —- iwé 5 =) (34) 
O O 


By uSing this expression in equation (9), integrating from 


weeeeeeo & — 1, (the hull/fin intersection) and adding the 


h 
contributions of the fins and thruster-rotor combinations, 
the complete turbulence forcing functions for the airship 
@em be found. That is, 


er - (Yg) , (55) 


tro 


Ih 
= + 
Xg | q(Yo) py 
0 J 
Likewise, the yawing moment about the nose is 


IL 
h a 
— ~d ; ‘| 
+) Bd(Yg), + 2 LO) 5 Yq) p, eee 


Nonose “ 
0 J 


v6 (ie, 


aud ig... “Sndse. ~ tom g 


and rolling moment is given by 


ee (37) 


a8 





and the non-dimensional equations can be expressed as 


oy = Gy Gy = Gy aeons (59) 

where y = Texp(iwt) = Texp (ikt) (60) 
eee 
2U, 


Gy, G 1, and G are the turbulent forcing functions for the 
YY ii ay 


vehicle. 
feeeerion Response Transfer Functions 
Using the functions of equation (59) as forcing 
E@metions on the right-hand Side of equation (53), and 


dividing through by y gives 


foe ee eye Cae | 
r 
Ra Cn. 
p/T | = | Si (61) 
24 
| O/T | Oo | 
byrr | fl, o 


Solution of this expression for specific reduced 
frequencies (kK), or spectral wave numbers (), and fixed 
Stable control gains (kK), allows calculation of the ex- 
pressions necessary for the solution of distributed force 


loadings and moments, by means of the following expressions. 
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ee. nee oo. £ ike 
Vig rt expi( lke), aT - exp (1kt) 
O O 
ee | 2B expres), » = 2 exp (ikt) (62) 
Uy r T 
Y . m~ 
y = r exp (1kt) 


feeuoOAD RESPONSE TRANSFER FUNCTIONS 
Once the motion response of the airship is known, the 
load response transfer functions can be calculated. 
ime Turbulence Loading 
These may be obtained by substituting equation 54 
into equations (9) through (13) and dividing by y (equation 


(60)). This gives, for example, for equation (9): 


Zz 
(aY_) U 
ac a ev exp (-in cosa_)dé (63) 
P ie 
where 2 = = eb 
U = 
O e 


Table II gives the complete list of load response transfer 
PiMmecLons . 
2. Motion Response Loading 
The aerodynamic-reaction loading may be obtained by 
replacing the motion variables, S|, eB etc., in equa- 
tions (14) through (20) with the Gieteencntdltinn motion- 


response transfer functions in equations (62). For example, 


equation (14) becomes 
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wh 2 pe APE 2 Gy gf) 
r > Giemane a Cit i 
Ae 2. 8 7 OFT ayia 
+ ou” a{ial¥ 2 2 ag-8) |, t= & Ky, pdé 
Cc Cc 
(64) 


The Inertial-Reaction, and Bouyancy loading transfer func- 


tions are similarly obtained from equations (21) through 





(27). For example, equation (21) becomes 
ST 20," R 2TV DR | 
—— ao roe = T its 120, im. 0 = Tem ®) | ean (6) 
t = : 


and equation (27) becomes 


(cies) 
=== [pgAd—é - gdm] cosa 2 ioe) 
r | 


Bie 


Saeammonrce, Bending and Twisting Moment Transfer 
Funct LONS 
The shear force loading transfer function is ob- 


tained using equation (29): 





l 
ayes) 
i rT ‘=] 
0 » 

where: - 

ay _ (dyY yn } (AY), i (eae. . (dY oy 

r r r i 
and 
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Ya) Os, 
x aa Ie 
ae i z J + ah] 
lr r r 


The bending-moment transfer function comes from equation (30) 








1 i (Yes) (Ya 

(a) (dN) wi. 

—-— = eee) ee eee. ; (1-19, || 3» 
0 iF 0 1 j=1 ae if 

(68) 


and finally the twisting-moment transfer function, from 


equation (31), is: 





JL i 
(dL_) (dL__) 
mil) _ -| h(E) (av). + | mh, | mg’ h 
r cm h A - 
0 
- (Yop (Yoo 
- } ho. a rl (69) 
j=1 J f r 


An important check on the analytical model is ob- 
tained by ensuring the net Shear Force, and Bending and 


Twisting Moments equal zero. 


SL) (om) (3) (Yo) (ay) 
in i CSP. ws F m s rs Ca 0 (70) 


r i: i P 














The moments are evaluated to the empennage assembly's mass 











center, Alon s! so that 
, 
BM (1 yo} ( You) (ae) 
om) s ee a gS 4, wees) f Ges 
ems S - 
Tr I i r [ 
(N_) (N_) 
ie Sipe (71) 
r Tr 





























TM(1__) (L,.) (L_) (Y_) (en) (Y_) o| 
cm's | ws ms , eas gas , ws. Cosi = 0 
r r r r r eo 
(72) 
These equations may be non-dimensionalized as follows: 
oe To ie) 
r ou” S 
Cam) 2 ap aa) 
= B r (74) 
T eU7se ~ 
O 
Coy (2) = 2 TM (1) 
~ —3o= ae) 
ig ou, Sc 


D. RESPONSE TO ATMOSPHERIC TURBULENCE 

Once the force and moment transfer function coefficients 
are known, the turbulence statistics can be applied to ob- 
tain estimates of airship lifetime and failure probability. 
When dealing with the lateral aerodynamic case, two different 
spectra must be considered--©,, and P55: This analysis will 
use the von Karman spectra as given by equations (7) and (6) 
respectively. As explained in Chapter II, the first is used 
when the flight direction is perpendicular to the mean wind, 
and the second when the direction is parallel. 

1. Root=Mean-Square Responses 

The root-mean-square response to turbulence of any 


system parameter can be obtained using its transfer function 
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multiplied by the spectrum (provided, of course, the spec- 
trum is Gaussian). For example, the rms shear force coef- 
frcient is 

co 


2 LZ 

s’ rms 2 iC. bao 

————————— = re —S dN (76) 
T 


NO 


2 
O 


QO ad 


Response to various conditions can be evaluated using appro- 
priate values for 9 and L in the equations for Daas and 
choosing either 211 ong 255 according to desired flight 
direction. 
2. Mission Analysis Method 

The mission analysis method 1s a technique for esti- 
mating a flight vehicle's probable lifetime. The method is 
based on the probability distribution of encountering tur- 
bulence in representative flight operations [Refs. 7 and 15]. 
It 1s assumed the total flight is a sum of Gaussian patches 
[Ref. 2]. The formula for calculating the number of ex- 


ceedences 1S 














-|x-x,o¢| -|x-x, 06!) 
N(x) = )tN_|p,exp nee + pi exp oe (a7) 
oy = Z = 
eee. biA 
il 2 
where: x = maximum structural value of bending, moment 


coefficient at a given station, 
Xof = value of x in one-g level flight, 
N(x) = average number of exceedences of the indicated 


value of x per unit time 


a4 





N = number of zero crossings of x per unit time 


A = ee ene oO Nee 


e = fraction of time in mission segment 


Pir Po = probability values from Table I 














bys b. = intensity levels from Table I 
also 
i (BM) a 
N = + MS 
Oo a (BM) as | 
where 
ou, “se 2 Com 
(BM) = |——- ——'1 9. .dQ (78) 
ms 5 - Beat 
0 
ou,*se ‘ oo 
(BM) = |-—_- —— 120. .d ress 
ms 5 - ie 
0 
-l 


The probable lifetime is then [N(x)] 

Shear force and twisting moment are analyzed in the 
same fashion using the appropriate transfer functions. 

3. Other Methods 

BeGeurizer and Hui [{Reft. 5] also include Failure- 
Probability analysis [{Ref. 3] and "Mil-Spec Storm" analysis 
in their paper. There are others in existence, such as the 
Design Envelope Analysis [{Ref. 7] that have been used for 
HTA flight, and are well documented. Because of this, they 


will not be included here, as techniques for using them are 
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the same once the transfer functions of response are known. 
The exact method used is up to the designer, based on his 


needs. 
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IV. NUMERICAL EXAMPLE 


In order to illustrate the lateral aerodynamic case 
developed in Chapter III, an example using the USS AKRON 
(ZR-4) is presented. The flight conditions chosen are: 


U, =o S fey) sec 


Alt = 1000 ft 
The velocity represents the maximum for the vehicle, and the 
altitude is typical of its operational range. In addition, 
a condition of neutral buoyancy (B-mg = 0.0) was selected. 
The geometry was taken from Freeman [Ref. 19] and the weight 
distribution from Woodward [Ref. 20]. With this information 
available, the inertial properties of the AKRON could be 
calculated using the method of Scholaert and DeLaurier [Ref. 
21] (see Appendix). The values obtained are shown in Table 
Tit. The Hull cross-flow and stabilizer efficiency factors 
are calculated using the method given in reference 6, and 
are: 


K Or ga22> 


n On Zale 


Ss 

The apparent-mass coefficients are from Munk [Ref. 22]. 
The stability derivatives are taken from DeLaurier and 
Schenck [{Ref. 18], and shown in Table IV. With these, the 


control gain, and the inertial and geometrical properties, 
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Seueteton 53 can ber solved to find the stable reots. This 
was done in reference 18 and the results are shown in figure 
8. Mode 4 (indicated in the figure) is characterized by 
roll, yaw, and sideslip of equal magnitudes not unlike the 
dutch roll mode of a fixed wing aircraft. Mode 5 is one of 
equal and opposite 8 and w motions with small ¢ perturba- 
tions. Mode 6 is a relatively high-frequency rolling motion, 
little affected by control gain. Modes 1, 2, and 3 refer to 
longitudinal aerodynamic modes [Ref. 18]. From this analy- 
sis, a control gain of 0.2 was found to provide the minimum 
seen le condition. 

The forcing functions were next obtained using equations 
54 through 60, and are plotted in figure 9. The peaks in 
all the curves occur at a wave number of about .008. This 
corresponds to a wavelength equal to the length of the air- 
ship. |S, | is significantly smaller than the others, as 
expected, due to the smaller moment arm through which the 
Seeetorce works in producing roll. 

With the turbulence forcing functions, equation 61 was 
solved to obtain the motion response transfer functions. 
These are shown in figures 10 through 14. Control gains of 
0.2, 1.0, and 2.0 were used to illustrate the effect of its 
variation. The most significant feature of these responses 
is the peak at a wave number of .008. This corresponds, as 
in the forcing functions, to a condition where the spectral 


component wavelength exactly equals the airship length. This 
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is the result predicted by Calligeros and McDavitt [Ref. 4] 
for the longitudinal case. DeLaurier and Hui [Ref. 5] also 
obtained this result, but only for cases of higher control 
gain. For the lateral case, the control gain does not sig- 
nificantly change this peak, although, for yaw and yaw rate 
(W and y), and to a lesser extent roll (%), the response at 
lower wave numbers is reduced. 

Finally, the load response transfer functions were cal- 
culated. The results are shown in figures 15 through 17 for 
a wave number of .009, and a control gain of 0.2. Complete 
results are given in the appendix. For the most part, the 
results yield no surprises. The magnitudes of the load 
response follow the general trend of the combined motion 
responses, thus the peak loads occur at wave numbers near 
.008. The location along the axis of the peak load varies 
as the magnitude of the motions increases, shifting aft in 
the case of shear and twisting moment, and to the center for 
bending moment. Again, the lack of significant change with 


Semeeol gain 1S apparent. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The analysis presented is an extension of the work by 
DeLaurier and Hui [Ref. 5], and is subject to the same re- 
strictions. That is, it is limited to small perturbations 
in order to allow a linear analysis usable with power spec- 
tral methods, and its ability to make precise poredictions of 
the loading when used with one of the methods that accounts 
for severe turbulence is questionable. Nonetheless, it is 
a valuable tool in understanding the response to an initial 
disturbance, and when employed as the aerodynamic input to 
tne various statistical methods discussed earlier, it can 
yield important design and operational insight. 

The limited effectiveness of the simple control model 
employed, which is typical of someone cuing his response to 
a compass, was demonstrated. It is suggested that an exami- 
nation of the effects of roll control and yaw rate feedback 
be made. This would allow a decision as to the feasibility 
of using control to provide gust alleviation. As discussed 
by DeLaurier and Hui, control gain made a large difference 
in the expected lifetime of an aie Si when considering only 
longitudinal aerodynamics. Undoubtedly, for the lateral 
case, even the yaw control used in this analysis will con- 
tribute to increased survivability due to the reduction of 


loads at low wave numbers. Thus, the next step is to employ 
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the statistical methods to discover how much change is 
realized. 

The case of combined longitudinal and lateral motion 
needs to be studied. No aircraft ever built has ever 
managed to fly through turbulence that is strictly one 
dimensional, as is assumed for this analysis. Coupling the 
two cases would give a much better idea of the true action 
of an airship in turbulence. 

Finally, some means must be found to establish the vera- 
city of this model, as well as that for the longitudinal 
case. To the author's knowledge, no investigation of the 
actual response of an airship to conditions of known turbu- 
lence has ever been made. This is, of course, a difficult 
project, considering the limited number of airships cur- 
rently available if full scale tests are to be carried out. 
Wind tunnel investigations, made in the various oscillating 
flow tunnels available, would be helpful. Until some tests 
are done, however, this tyve of analysis must be considered 
only for its qualitative aspects as opposed to its quantita- 


tive predictions. 
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TABLE III 
GEOMETRICAL AND INERTIAL PROPERTIES 


OF THE USS AKRON (ZR-4) 


total volume = 7,382,400 ft 


eD..0) fic 


ey one 


Cc 


i 364.24 £t 
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eo OO te 


ly =—ooo.0 lL £t 


mass = 17,039 slugs 


buoyancy = 548,642 lb 


. Le 
I, = 38,685,100 slug-ft 
I, = 471,799,000 slug-ft? 

_ 2 
I, = 102,129,000 slug-£t 
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TABLE IV 


Sloane, DERDVALTLVins OF THE USS AKRON (ZR-4) 


neutral buoyancy, we = 123 ft/sec, ALT = 1000 ft 
Cy, -0.7224 Cys -0.9863 
Cy -0.3418 Aya 0.0586 
Cys ~0.0648 Cys =Q0o 13 
Cne rer ng Cng O. 9293 
Ch. = -0.2352 Che = -0.0991 
Cn, = -0.0150 Cns 0.0028 
Cl, = -0.0322 Cig = -0.0456 
Cie -0.0153 Cle O00 3 
Cy, = -0.0066 Cy, = -0.0042 
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Figure 1. The (1-Cosine) Gust Shape 
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Figure 3. Elementary Spectral Components 
in Two Dimensions 
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Figure 4. Typical Power Spectra of Vertical Gust Velocity 
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Figure 5. Measured and Fitted von Karman Spectra 
Of Vertical Gust Velocity from Severe Storm 
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Figure 7. Schematic of Buoyancy Forces and Moments 
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USS AKRON (ZR-4) 
Neutral Net Buoyancy MODE 4 eT 
US = 123 ft/sec 
wie. — L000 £t. 
0.8 
k = 1.0 
C 
0.6 
0.4 
Sm a 
cla = MODE 5 ee 
/ / / k, = 1.0 
alee 1.0 -0.8 -0.6 -0.4 “0.2 0. 
Cr 


Figure 8. Lateral Root-Locus of the USS AKRON (ZR-4) 
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